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Abstract 

The problem of interest in this article are waves on a layer of finite depth 
governed by the Euler equations in the presence of gravity, surface ten¬ 
sion, and vertical electric fields. Perturbation theory is used to identify 
canonical scalings and to derive a Kadomtsev-Petviashvili equation with 
an additional non-local term arising in interfacial electrohydrodynamics. 
When the Bond number is equal to 1/3, dispersion disappears and shock 
waves could potentially form. In the additional limit of vanishing elec¬ 
tric fields, a new evolution equation is obtained which contains third- and 
fifth-order dispersion as well as a non-local electric field term. 


1 Introduction 

Classical water wave models leading to equations such as the Kadomtsev-Petviashvili 
(KP) equation have long since been used to understand the nonlinear phenom¬ 
ena and interactions of waves and are useful in the theoretical basis for further 
studies. The main interest here is to extend the results in [8] to the three dimen¬ 
sional case. The importance of interfacial electrohydrodynamics phenomena has 
been highlighted in many different cases. 


2 Set-Up and Governing Equations 


Consider a perfectly conducting, inviscid, irrotational and incompressible fluid 
(region 1) bounded below by a wall electrode aX y = —h and bounded above 
by a free surface y = r](t,x^y), here h is the mean depth of the surface. The 
fluid motion is described by a velocity potential (p{t, x, y, z) satisfying Laplace’s 
equation in region 1. Surface tension with coefficient cr and gravity, g, are 
included. The region y > ri(t,x,y), denoted by region 2, is occupied by a 
hydrodynamically passive dielectric having permittivity Cp. It is assumed that 
there are no free charges or currents in region 2 and therefore the electric field 
can be represented as a gradient of a potential function, E = — W. A vertical 
electrical field is imposed by requiring that V ^ —E^y as y —>■ oo, where Eq 
is constant. The voltage potential satisfies the Laplace equation. On the free 
surface the Bernoulli equation holds: 
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The pressure p is obtained through the Young-Laplace equation: 


[n • T • n]2 = crV • n 


Where the stress tensor is given by: 


T = - 


pSij + Cp ^EiEj 


The unit normal is given by: 

~ ^ [dxV,dy-q,-l) 
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The governing equations are then: 
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V ^ —Eqz as z —>■ oo 
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3 Linear Theory 


We move into a reference frame which is moving with the fluid, 
the time derivatives from the governing equations. 

This will remove 


3.1 Finite Depth 

The scaling we use for the finite depth case is: 

X = hx, y = hy, z = hz, p = hr) 

V = EohV, P=^p 
h 
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There are two dimensionless parameters which come out of the non-dimensionlisation 
process: 


B = 


ph^g 


Eb = 


eE^ph 


(13) 


a a 

The parameters are called the Bond number and electric Bond number respec¬ 
tively. The variables are expanded as: 


p = epi -I- o(e) 
r\ = erji + o(e) 

(fi = Ux + eipi+o{e) 
V = —z + sVi + o(e) 


The variables Vi,T]i,(pi,pi are written as inverse Fourier transforms: 

1 


Vi = 


Pi = 


m = 


(2^)2 

1 

1 


E'^^^+^^yA{k,l,z)dkdl 

e^^i^+ayB[k,l,z)dkdl 

E’^‘^+^iyC{k,l)dkdl 


The expression for the free surface is then given by: 


, , 1 /* 

= ^ k'^U^ - {Bp - Ebp^ + ’ 

where p = y/k"^ + B. The stability relation is given as: 


(14) 
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Typical free surface profiles predicted by (1^ for Fit, = 1.5, i? = 2 and U = 0.5 
in figure (13.11) . 


3.2 Infinite Depth 

The difference between the hnite depth and the infinite depth is the boundary 
conditions, these new boundary conditions are given by: 


Using the expansion: 


dp 

dx 

—>■ U as z oo 
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Figure 1: Linear Waves Profiles - Finite Depth 


The same approach can be done for the infinite case, by writing the variables 
as inverse Fourier transforms, it is possible to write the free surface as: 


Vi = 


WP 






Akx+lv)dkdl 


(29) 


Where p is the Fourier transform of the pressure. Equation (1291) can be made 
clearer by the introduction of the following variables: 


a = U^g /3 = U'^g ^1, x = xgU y = ygU = gpU ^ (30) 


To obtain: 


where: 


f)i{.x,y) = 


1 f 

— v + piA — ^21^^ 


dad/3 


Pi = 




pU^ pU* 

Figure ([5]) show the waves for infinite depth for values pi = 1 and p 2 = 2. 
When /3 = 0, the denominator in the integrand can be arranged to be a perfect 
square provided that pi and p 2 satisfy a certain relation 




(31) 


(32) 


1 + /i2 — 2^/p2. 


(33) 


4 Weakly Nonlinear Theory 

Previous work in this area has been carried out by Katsis and Alylas mi 
whereby a two fluid scenario with interface was given by z = x, y). Analysis 
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Figure 2: Profile of 2D Infinite Depth Wave 


showed that the resulting equation was the same as that obtained in this section, 
despite there being a different method and different coefficients. For the weakly 
nonlinear theory, use the following scaling ([I],[H],IZ]): 


A ~ ^ gXa ^ 

X = Xx, y = yy, t = — t, rj = arj ip = - ip 

Co Co 

V = XEoV = hz = Ai, 

Define three parameters mm-- 


a = 


h' 


^ A2 ’ ^ y? 


The governing equations become: 

1 d'^p d^ip 


d^V d^V d^V „ . A 
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dlfjil + a^P^{dyf]Y) + 79 ? 77(1 + a^P{d±fjY) - 2a^p-^dxf]dyfjd±dyf] 


(1 + a^P{d±f]Y + a^P^{dyfjY)‘i 


Where: 


B = 


pgK^ 


Eb = 


pgh 


(43) 

(44) 


The non-dimensional parameter B is the Bond number. Note that Et can be 
thought of as the ratio of two things, CpE^/p and gh. The quantity gh has the 
units of speed^ and therefore so does CpE^/p. The speed is therefore 

characteristic to the system of interest, the quantity \/Eb can be thought of as 
the electric Foude number. The constant was calculated by using the solution 
ip = fj = 0 and V = —z to show that C = —Eb/2a. To begin with restrict the 
attention to the classical shallow water scaling by taking: 


a = /3 = 7 = e<Cl, T = ei, X = x — i (45) 


From here on in, drop the hats. The asymptotic expansions for the variables 
are: 


ip = ipQ+ Eipi + e^ip2 + o(£^) (46) 

g = gQ + egi+ o(e) (47) 

V = —z + e^Vi + o(^e^^ (48) 

P = epi (49) 


The expansion for ([35]) comes from examining Equation ([57|) can be solved 
as: 

{z + lfd'^PQ ^ 2 r + 1)"^ ^^¥^0 {z + l^d’^ifo 

^ 2! dX^^^ [ 4! dX^ 2! dy"^ 

Where the equation shows that ipo = ipQ{T,X). The electric term in the 
Bernoulli equation is: 


+ o{e^) (50) 


^ r 1 ^ dVi 

T 


(51) 


The hrst part of this equation cancels the Bernoulli constant and then this 
leaves a term of order £ 2 , so to include this into the order e equation, the 
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electric Bond number is scaled according to Ei, = Eby/s- This makes the 0(e) 
Bernoulli equation become: 

?7i = Bd\ - dripo + -pi+ Ebd^V - (52) 

Moving on to the free surface equation (HTl) . the 0(1) equation gives rjo = dxfo, 
which was known previously and the 0(e) equation is: 

^dxPo - - ]^'nodxPo = drm - dxm + dxfodxVo (53) 

Inserting the expression for rji coming from (1521) shows that: 


^djcipo-dlpo = 2dTdx<fo+^ Q 
In terms of rjQ the equation is: 


9xPo+dxPiE—dxij]Q)+EbdxdzV — 0 

(54) 


d dr]o 1 I'1 \ d^r]Q dpi 3 dr]o Eb d'^V 

dX dT ^ 2\3 ) dX^ ^ dX ^ 2 '^° dX ^ 2 dXdZ 

The 0(1) equation for (I!?R1) needs to be solved, it is given by: 


1 d'^ijQ 

2 dy'^ 


= 0 


(55) 


dX^ ^ dZ^ 


(56) 


This equation requires a boundary condition in order to write down a solution. 
All boundary conditions are zero except the one at Z = 0. The boundary 
condition can be found by expanding (IdHll to yield: 


dVi dp 


(57) 


Integrating this equation and assuming that the free surface dies off at infinity 
shows that Vi = 77 on Z = 0. Equation (l56|) is solved by use of a Green’s 
function. The Green’s function for a 2D Laplace equation in the upper half 
plane is: 

{X - X'f + [Z - Z'f ^ 

{X - X'Y+ [z + z'y\ ^ ’ 

and use of Green’s second identity: 


g(X,ZlX',Z') = —log 


[ gX^Vi-ViX^gdT= [ g{n ■V)Vi - ■V)gdX 
JD JdD 


(59) 


Inserting the boundary condition into (1591) shows that the solution is given by: 



TT 


Zpo 

{X - A')2 + Z2 


dX' 


Then the required result is: 


dVi 

SI “ ( m ) 


(60) 


(61) 
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So inserting (ICTl) into (1551) shows that: 
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dr]Q 

dT 
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2 V3 


- B 


d^Vo , dpi , 3 dr]o , h ^ (d'^r]o 
dX^ dX 2"^° dX 2 I dX^ 


15% 
2 dy"^ 


= 0 


In terms of dimensional variables the equation is: 


(62) 


dx Co dx 2h^dx 2 \3 / dx^ pg dx 2 \ dx"^) 



(63) 

When p = Eb = 0, (1551) reduces to the standard KP equation. In the 2D case, it 
has been shown that if i? < 1/3, then “generalised” solitary waves are possible. 
The equation which we have derived is exactly the same as equation (2.7) in 
[13] but with Eb replaced with the ratio of densities of the fluids. So there 
appears to be a link with interfacial flows, at least on the level of equations 
of motion. There is different behaviour around B = 1/S and the next section 
details the derivation of a new equation in this case. Note that a complete 
analytical derivation of all terms in the equation has been given whereas they 
simply add the nonlinear term in m- The dispersion relation can be easily 
derived for the problem, it is given by: 


uj = p\g 


cEq O' 2 \ 

-/i H —p tanhh^ 


(64) 


where p = and k and / are the wavenumbers. Equation (1641) can be 

expanded in powers of k and I to obtain the linear portion of (1631) . 


5 Analysis around 5 = 1/3 

In this section, the scaling in (1551) is replaced by: 


a = e^, /3 = e, 7 = £^, T = X = x — i 

The expansions are now given by: 

(65) 


= ipo + eipi-\-e^ip2++ o{e^) 

( 66 ) 

V 

= 7o + £Vi + + o(e 2 ) 

(67) 

V 

= — Z + Vi + 0 

( 68 ) 

p 

= e^pi + o{e^) 

(69) 

B 

= — + eBi + o{e) 

(70) 


The governing equation for ip is solved in exactly the same fashion as in the 
previous case, the solution is given by: 


P = <fo{T,X) 


(Z + 1)2 d\o 

^ 2 ! dX^ 
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, 2 {z + 1 )^ d^ifo 

4! dX* 


(71) 












Only the derivative of the third term needs to be calculated, this is: 


dipz _ 1 _ d'^ipo 

dz ~ 5! dX^ dy^ 


(72) 


The 0(1) Bernoulli equation once again shows that dxfo 
part of ((1^ is: 

1 d^Vo _ 

3 dX^ ~ dX 


r]o, the 0(e) 
(73) 


The term ipi has been calculated in (ED, the O(e^) part of the free surface 
equation m is: 


dryo — 9x^2 + dxifiodxVo = —iJodx^PQ + dzips (74) 


To find ? 72 , the 0(6:^) Bernoulli equation is used. The same rescaling as 
before is used with the electric Bond number, Eb as Eb = EbS'^ in order to 
include the electric term in. The equation becomes: 

BidxVo + \^xVi = dripo - dx<f2 + ?72 +Pi + ]^{dx<^oY - EbdzVi (75) 
Inserting this into (1741) shows that: 


29 t ? 7 o - Bidx-rjo + ^rjodxVo + '^^xVo+ 

+ dyipo + dxpi + EbdxdzVi = 0 


(76) 


In the previous section, the term dxdyVi was examined before and the result is 
just quoted, dxdyVi = J^{dxr]o), making the final equation: 
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dX 


dyo 1 d^r]Q 
dT 90 dX^ 


3 ^770 Bi d^po 


2^°'dX 




2 dX^ 

1 d^r]o 

2 dy"^ 
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In dimensional variables the equation becomes: 


A 

dx 


1 dr] dr] d^r] 3 dy f 

Co dt dx~^ 90 dx^ 2h^dx 2 \ 3/ dx^ 


{ Tr4\ + +xzr4=0 






2pg dx 

2 dy"^ 


(78) 


6 Fully Nonlinear Results in Infinite Depth 

Only the infinite depth case in considered for simplicity, but an extension to 
include the finite depth is possible. The conditions at infinity are 

V ~ —EqZ, for z —oo. 

d<l) , r 

— ->■ I, tor z —>• —oo. 

ax 
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The interested is in steady waves travelling with a constant speed U, so choose a 
frame moving with the wave and non-dimensionalise the equations by using the 

unit velocity U and the unit length L = -Also V is non-dimensionalised 

pU‘^ 

using EqL. The non-dimensional parameters are 


^ ^ _ 

pu^ “ pm 


These are the exact same parameters obtained in the linear case with pi = Eb 
and p 2 = This is very useful as it can be used to determine the rage of 
validity of the linear solution. 

Set = 0 at the free surface z = r]{x, y) and introduce V = V + VqZ, which 
satisfies y —0 as z —oo. The Bernoulli’s equation becomes in this steady 
frame 


1 

2 





——+ar]—K— 


2 



a/i+^ 





(79) 

where n is the downwards unit normal (see equation (4)), and K is the curvature 
given by 


Vx 


Vy 

^J^ + vl+vl_ 

X 

^J^ + Vx+Vy_ 


By applying the second Greens identity for (j){x,y,C,{x,y) — x in the region 
z < ? 7 (x, y) and for V in the region z > r]{x, y), we obtain the boundary integral 
equations 


(80) 

= -ZZ^,(xp,„,) a . V dSp + ZZ^ 

(81) 

where P and Q are points on the free-surface Sp with coordinates {xp, yp, C{xp,yp)) 
and {xQ,yQX{xQ,yQ)) and 


RpQ = \PQ\ = \J {xp - xqY + {yp - yq? + (C(a;p,2/p) - C,{xQ,yQ)Y. 


The equations (I29|)-(I8T|) are desingularized and projected on the Oxy plane, 
and then integrated numerically (see Landweber & Macagno 1969, Forbes 1989, 
Parau & Vanden-Broeck 2002 for details). For most of the computations 40 
grid points are chosen in x direction and 40 points in y direction, with a grid 
interval of dx = 0.8 and dy = 0.8. The fully-localised solitary waves near the 
minimum of the dispersion relation have been computed which corresponds to 
the region 

/ p \ 2 


a > 


\ + Eb 
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It should be mentioned that for Eh = 0 this minimum corresponds to a = 1/4, 
which is the critical point for gravity-capillary waves (see Parau et al. 2005). 

A typical profile is given in figure [31 



Figure 3: Fully Nonlinear Infinite Depth 


7 Discussion and Conclusions 

In this paper an expression for the free surface of a perfectly conducting fluid 
under the influence of a vertical electrical held in the linear and weakly nonlinear 
cases has been derived in a simple, systematic way. It has been shown that, in 
general, to be an extension of the results for the simple Kadomtsev-Petviashvili 
(KP) equation and that this seems to be the generic expression for almost 2D 
equations for waves. A new equation has been derived (17811 which incorporates 
a number of different physical phenomena, surface tension, electric helds and 
when the Bond number is close to 1/3. Although the equations which have been 
derived here have appeared elsewhere ([T3],[I1]) with either different coefficients 
(as in m) or with some terms absent (as in M) a simple derivation of the 
results has been obtained whereas in other works special assumptions have to 
be made with the coefficients ([2]) whereas the coefficients come out naturally 
in the analysis presented here. 
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